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Abstract 

We theoretically explore the notion of nonreciprocal near-zone manipulation of electromagnetic 
fields within subwavelength plasmonic nanostructures embedded in magneto-optical materials. We 
derive an analytical model predicting a strong, magneto-optically induced time-reversal symmetry 
breaking of localized plasmonic resonances in topologically symmetric structures. Our numerical 
simulations of plasmon excitations reveal a considerable near-zone power flow rotation within such 
hybrid nanostructures, demonstrating nanoscale nonreciprocity. This can be considered as another 
mechanism for tuning plasmonic phenomena at the nanoscale. 
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Recent advances in plasmonics - nanoscale electrodynamics of metals - have provided 
roadmaps for significant miniaturization of optical devices and components. Consequently, 
electromagnetic fields can be confined, enhanced, and manipulated in length scales as low as 
few nanometers around plasmonic structures jl|. Such a dramatic field confinement led to 
a number of breakthroughs reported lately, including nano-lasers j2|, medical treatment jsl, 
subwavelength nano- antennas 0,3], observation of novel magnetic resonances 6, 7| and some 
other applications jsj. However, further scientific development and technological expansion 
in the fields of plasmonic optics and nanophotonics may in general be accelerated by the 
ability to tune actively the optical response of plasmonic strutures and their field distribu- 
tions at the nanoscale. Several techniques developed in recent years such as active all-optical 



control 
ifications 



Q,U,li8l 



nonlinear self-tuning 



12 



13|, waveform shaping 14|-|l6| and structural mod- 



are regarded as possible pathways for the nanoscale field manipulation. 



In particular, in Refs. 0, ll8| a structural design and arrangement of plasmonic oligomers was 
exploited for precise control over the resonances and corresponding near field distributions. 

Mixing magneto-optical (MO) materials with plasmonic structures may provide another 
mechanism, and an additional degree of freedom, in tailoring the light-matter interaction in 
the vicinity of plasmonic structures. It is well known that the MO materials may possess 
nonreciprocal response as a result of breaking the time-reversal symmetry in optical phenom- 



ena 



19|-|2l|, - a property that may be exploited for signal handling and manipulation, and 



for photonic circuits 2ll-l24j. However, the magneto-optical response of naturally occurring 
materials is usually very weak 25|, and may not be adequate for some potential applica- 
tions. Recently an incorporation of MO materials with plasmonic structures has gained 
a significant attention, due to possible enhancement of magneto-optical activity in highly 



localized fields |26l . 



27|. Enhancement of such macroscopic (i.e. far field) nonr ecip rocal 



effects as Kerr and Faraday rotations 



26 



30l |. and transverse Kerr effect 27|, |3l| have 



been demonstrated. However, to the best of our knowledge, the microscopic (i.e. near-field) 
nonreciprocal optical response and its enhancement have not received as much attention. 

In this Letter, we theoretically explore and analyze the concept of nonreciprocal ma- 
nipulation of near-zone optical field and nanoscale power flux by mixing magneto-optical 
materials with plasmonic nanostructures, and demonstrate, using numerical simulations, 
a significant enhancement of the nanoscale nonreciprocal response. In particular, we study 
field distributions at plasmonic resonances in such structures and show that magneto-optical 
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FIG. 1: (color online) Schematic of the two structures studied in the Letter: a) A plasmonic 
nanorod in an equilateral triangular magneto-optical cavity surrounded by a free-space; and b) a 
collection of nanorods embedded in an unbounded MO medium. Magnetization field shown by a 
red arrow is directed along the main cylinder axis. 

activity may lead to a strong coupling between the degenerate eigenstates corresponding to 
the same resonant frequency, resulting in the formation of highly localized rotating eigen- 
states. We analyze the plane wave excitation of corresponding solutions that depend on the 
strength of magneto-optical activity, and reveal a pronounced power flux circulation. We 
discuss typical structures where the boosting of near-field power flow circulation is observed, 
and give a brief insight to the possible experimental realization and potential applications 
of such observed effect. 

We begin our study with the a generalized two-dimensional (2D) eigenfrequency analysis 
for an arbitrary 2D plasmonic (i.e. metallic) structure embedded in a magneto-optical 
material. Fig. [1] shows geometries of the problem with the z axis of a Cartesian coordinate 
system being parallel with the axis of geometries. The magnetization of the MO material 
may be parallel or anti-parallel with the z axis (i.e. the Voigt configuration). We search for 
TM electromagnetic field distributions in the (x — y) plane, for which the electric field has 
X and y components and magnetic field only z component. In this case the MO material 

3 



response is given by an antisymmetric relative permittivity tensor 25|: 
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where Smo and e_[_ are diagonal components of dielectric permittivity, a is the off-diagonal 
component of the permittivity tensor responsible for the "strength" of magneto-optical ac- 
tivity of the media. Typically a at optical frequencies is very small and usually is of the 
order of 10~^ or even smaller. The wave equation for the electric field in the (x — y) plane 
can be written in the following operator form: 
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where and Ey are the electric field components, u is the radian frequency, c is the free- 
space speed of light, e{x, y) = Smo in MO material and e{x, y) = Sm in metal, and ({x, y) = I 
in MO media and C{x,y) = outside of it. Obviously, last term in Eq.([2]) associated with 
the MO activity, leads to a nonreciprocal coupling between the E^ and Ey electric field 
components. When a = the Eq.(l2]) reduces to an ordinary wave equation. Solution of an 
eigenfrequency problem in a nonmagnetized case yields a set of resonant frequencies with 
corresponding set of eigenstates („ m)' -^y (n m))'^ where the parameters n and m are the 
modal parameters related to the radial and azimuthal variations of the mode, respectively. 
In general, the mode degeneracy may be present, implying that d {m = l..d) eigenmodes 
may possess the same given resonant frequency a;°. Note that, for a nondegenerate sys- 
tem the magneto-optical activity causes the interaction between eigenmodes corresponding 
to different resonant frequencies u^, however, when these resonances are well pronounced 
and well separated from each other such intermode interaction can be considered negligi- 
ble. At the same time the maximum magneto-optical interaction is expected between the 
degenerate states corresponding to the same resonance, and in further analysis we focus on 
this case only. Considering that a << Smo, we apply the perturbation method for general 
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case of d times mode degeneracy at a a given resonance and study the interaction be- 
tween these degenrate modes. We consider the series expansion of the eigenfrequency and 
corresponding eigenmodes in terms of powers of a as: (cOnY = i'^nY + '^i^nY + 
Eu,n = X] '^m-E'°(„m) + '^^u n + whcrc u stands for the coordinate x or y, and c^'s are 
unknown complex coefficients. Substituting these relations in Eq. and applying standard 
methods of perturbation theory, we obtain: 



+i 



f 1 «.)*«J + 








(3) 



where the summation is taken over the index m only (m = l..d), ipn,m — i-^x (n m)' (n m))^ i 
and u\ is the first order correction to the resonant frequency a;° due to the presence of MO 
activity, i.e. when a is non-zero. 

Assuming an idealized lossless system, we obtain real valued solutions if)^ ^. in the zero-th 
order approximation (i.e a = 0). In this case it is possible to show that the last term under 
the integral in Eq.Q vanishes for m = k, whereas the first term being the scalar product 
of the eigen-modes is zero for m ^ k due to the orthogonality. Taking these into account 
one can show that Eq.(l3]) reduces to a eignevalue problem for a d x d matrix with Cm being 
the eignevectors of this matrix and a;^ unknown eigenvalues. The solutions of such matrix 
are defined by its zero determinant. Eq.(l3]) shows that MO activity leads to the coupling 
between the degenerate eigenstates, which depends on the cross-product between the k and 
m states, i.e. the last term in Eq.([3]), and results in the formation of a novel set of complex 
eigenstates (^mi^n m corresponding to the l-th eigenvalue of Eq.Q. The physical meaning 
of such coupling is that the back-and-forth energy beating between the A;-th and m-th states 
is established. 

For a doubly degenerate case the determinant of Eq. ([3]) is written in a simple way: 
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FIG. 2: (color online) Intensity profiles of magnetic field for the degenerate states in the two 

structures at their second order resonance (n = 2) and cross-product between these states for a 

nanorod in the MO cavity (a,c,e) at ojq = 404THz, and the collection of nanorods embedded in the 

MO medium at ujq = 336THz (b,d,f), respectively. 
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where < • >= J •dS. Considering that < {i^niY >—< ('0n2)^ > '^^ derive the expression 
for the frequency sphtting: 



- " I' ^ i;;: 2 < + > ) 

here we have dropped indexes n and for simphcity. In this case the eigenvectors of Eq. (j3]) 
for d = 2 are complex conjugates, which means that the solutions of Eq.(l3]) to the fist order 
approximation are E± = {E^i ± iE^2, Eyi ± iEy2) with corresponding eigenfrequencies u+. 



In analogy with clockwise and counterclockwise plasmonic states of a single nanorod 32 1 



these new solutions E± can be classified as counter-rotating eigenstates. IN 
solutions were obtained previously for a photonic crystal circulator 23 



ote that similar 



2J|. When the 



frequency split |a;+ — a;_| is relatively large, these states can not be excited with the same 
strength by an incident electromagnetic field, and therefore the superpositions of these states 
gives rise to a near-zone field with a rotating near-field pattern. It is important to note here 
that for doubly degenerate structure possesing a 90 degree rotational symmetry between 
the eigenmodes (simplest example of such a structure would be a single plasmonic nanorod 
surrounded by MO) the cross product < {ExiEy2 — EyiEx2) vanishes, implying that the 
reported mechanism of mode coupling and corresponding symmetry breaking do not take 
place. 

The specific type of modal degeneracy depends on the structure design. Generally the 
degeneracy is observed for structures with some rotational symmetry. As two examples, 
we consider here a 2D plasmonic nanorod in a MO 2D cavity surrounded by air, and a 
collection of parallel nanorods immersed in a MO background (see Fig. [1]). In particular, we 
study a plasmonic nanorod with radius R = 50nm centered in a MO equilateral triangular 
cavity with side 450nm, surrounded by air. Fig. dJ^a), and a plasmonic nanorod with radius 
R = 50nm surrounded by four other nanorods with smaller radii R2 = lOnm symmetrically 
centered around the first nanorod, 80nm away from its axis. Fig. WO^)- Without loss of 
generality, here we assume that the lossless metal dielectric constant is fixed at Sm = —10, 
and we consider our MO material to be a Bismuth Iron Garnet (BIG) with Smo = 6.25 and 
a as a free parameter for our parametric study. 

First, we consider no MO activity (i.e. a = 0, for zero magnetization) and search for 
the eigenfrequencies and corresponding modal distributions using numerical simulation with 



7 




FIG. 3: (color online) (a) Frequency splitting |a;+ — vs magneto-optical activity parameter 
a and (b) net azimuthal power flow around the main cylinder, < >. Curves 1-3 correspond 
to the nanorod in the MO cavity, the collection of nanorods embedded in MO medium, and the 
single nanorod surrounded by uniform MO, respectively. Curves 4,5 the same as 1,2 but for a lossy 
nanorods. 



COMSOL Multiphysics, a commercially available finite-element-based simulation software. 
The eigenfrequency analysis shows that both structures possess two doubly degenerate modes 
at the second order plasmonic resonance [n = 2). The corresponding distribution of the 
magnetic field intensities for both degenerate states are shown in Figs. |2](a,c,b,d). One 
can see that the fields are highly localized in an area less than 100 x lOOnm^. The cross- 
products of these degenerate modes are shown in Figs. [2](e,f). Clearly, the cross-product 
is azimuthally symmetric and monotonically decaying along the radial direction, implying 
that overall integral < {ExiEy2 — EyiEx2) >\mo would be nonzero and therefore a strong 
magneto-optical coupling will exist between the corresponding modes, if the MO activity is 
introduced. When the structure is magnetized, i.e a 7^ 0, the two counter-rotating states 
with u± eigenfrequencies are formed, as was mentioned above. In Fig. |3](a) we plot the 
values of the frequency split — as a function of MO parameter a, using Eq.(j5]). We 
note that even for small values of a, the frequency split is significant and is up-to 10 times 
higher than that reported in Ref . 23| . Noticeably the coupling strength for both structures 
studied in this Letter is almost analogous, which can be attributed to the similarities in 
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the topology of near-field distributions, see Fig. [2l It is important to compare the reported 
mechanism of symmetry breaking with that of single plasmonic nanorod {R = 50nm) sur- 
rounded by magneto-optical media. Although the nanorod has two degenerate states at 
each plasmonic resonance {n = 1,2,...) the cross-product between them is zero, implying 
that the reported mechanism of symmetry breaking is not applicable for such nanorod so 
that the corresponding frequency splitting is two orders of magnitude smaller. Fig. 3 (a) (note 
that the MO activity still leads to the symmetry breaking between these degenerate modes 



associated with an interface effect 



33j). 



Next, we study the excitation of the predicted rotating states. In each of the geometries 
studied here, we consider a plane wave incident on the structure, as schematically shown in 
Fig. [TJ We numerically simulate this problem using COMSOL Multiphysics software. We 
trace the near-field distributions at these observed resonances and analyze their variation as 
a function of the MO parameter a. In Fig. |4] we present distributions of magnetic field and 
the Poynting vector for a nanorod in a triangular MO cavity at its second order plasmonic 
resonance (n = 2). We observe a strong distortion and one way-rotation of the near-field 
pattern and its power flux, when MO activity is present. The Poynting vector, describing 
local power flow, experiences a dramatic change in its behavior as one increases a. In 
particular, for an a = case we observe a mirror-symmetric power flow around the nanorod. 
Fig. m^a), i.e. there is no preferred direction of rotation for this power flow. However, when 
the MO parameter a is non zero, we note a strong power flow circulation around the nanorod 
as shown in Fig. ID^b). Similar behavior is found for the structure shown in Fig. [I](b), where 
again the magneto-optical activity leads to the breaking of mirror symmetry of power flow 
around the structure (not shown here). We note that the analysis of fleld distributions at the 
flrst order plasmonic resonances {n = 1) does not reveal any pronounced energy circulation, 
although both geometries are doubly degenerate at this resonance as well. The latter is due 
to the zero cross-product between the modes. 

In order to quantitatively determine the strength of such Poynting vector circulation 
we evaluate the net azimuthal power flow by inscribing the structure into a mathematical 

L 2n _ 

circular region and then we calculate < >= \J J (Pd(j)dr\, where is the azimuthal 



unit vector, P is the Poynting vector, and L = 150nm is the radius of our integration 
domain arbitrarily selected. In Fig. I^b), we plot this value for both of our structures. 
We observe almost linear increase in the net power flow circulation with the increase of 
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FIG. 4: (color online) Magnetic field (colormap) and Poynting vector (arrows) distributions for 
a plasmonic nanorod in a triangular MO cavity at the second order plasmonic resonance (u = 
404THz) when a = (a) and a = 0.1 (b). 

MO activity. Clearly, the rotational power flow as a function of a correlates well with the 
frequency split vs a shown in Fig. |3]^a). Owing to the nonreciprocal nature of the problem, 
the rotational direction and net power flow of circulation are dependent on the direction of 
the magnetization. We also study the plane wave excitation of a single nanorod surrounded 
by a magneto-optical medium, and analyze the corresponding net azimuthal power flow 
circulation, see Fig. [3t^b). As can be seen, the overall circulation in this case is almost two 
orders of magnitude smaller, ensuring that the near-fleld rotation is negligibly weak, which 
agrees well with frequency splitting predictions shown in Fig. |3t^a). 

Finally, we study the influence of losses on the system dynamics. In particular, we trace 
the power flux circulation in a lossy system taking into account the metal loss, i.e. consider- 
ing that lm{em) =0.1. Corresponding net azimuthal power flow is shown as dashed curves 
in Fig. [3]^b). Evidently, the nanocirculation of the power flux around the plasmonic nanos- 
tructures is well preserved in these lossy systems, and is one order of magnitude higher than 
that for a single lossless cylinder in unbounded MO medium. The predicted phenomenon 
of nanoscale circulation of power flux due to combination of the plasmonics and MO ac- 
tivity may have interesting impacts on various scenarios ranging from optical tweezing and 
nanoscale paticle manipulation and rotation, to nanoantennas, and nanocircuits. 

In conclusion, using analytical and numerical study we have shown that incorporating 
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magneto-optical materials with plasmonic structures leads to a significant enhancement of 
near-field nonreciprocal response. We have demonstrated that in the presence of MO activ- 
ity the modes excited at plasmonic resonances in certain types of geometries couple strongly 
with each other resulting in a formation of nonreciprocal rotating states. The plane wave 
excitation of two plasmonic geometries studied here has revealed a strong nearfield power 
flow circulation. This near-zone subwavelength power flux circulation and rotation may be 
exploited for design of next generation tunable nanoscale plasmonic devices. The principles 
of nanoscale optical field manipulation might pave the way for a wide variety of potential ap- 
plications in tunable sensing, active plasmonic elements, particle manipulation, and others. 
This work is supported in part by the US Air Force Office of Scientific Research (AFOSR) 
grant number FA9550-10-1-0408. 
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